Theory of tunneling spectroscopy for chiral topological superconductors 
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We study the charge conductance of an interface between a normal metal and a superconducting 
quantum anomalous Hall system, based on the recursive Green's function. The angle resolved 
conductance ^{ky,eV) with the momentum ky parallel to the interface and the bias voltage V 
shows a rich structure depending on the Chern number J\f of the system. We find that when the 
bias voltage is tuned to the energy dispersion of the edge mode, eV — i?edge(fcy), the angle resolved 
conductance 7(fcy, i5edge(fci/)) shows a pronounced even-odd effect; the conductance vanishes for 
J\f — or 2 while it takes a universal value 2e^ jh for M —\. In particular, in A/" = 2 phase, we find 
that the conductance ^{ky, Eedge{ky)) becomes zero due to interference of two degenerate Majorana 
edge modes, although the corresponding surface spectral weight remains non-zero. 

PACS numbers: 74.45.-l-c, 74.50.-|-r, 74.20.Rp 
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I. INTRODUCTION 

It is well known that Andreev bound states (ABSs) 
are generated at the edge of unconventional superconduc- 
tors where the pair potentials change sign on their Fermi 
surfaces.—"— Nowadays, the ABSs in unconventional su- 
perconductors have been recognized as important topo- 
logical objects in condensed matter physics. Supercon- 
ductors hosting topologically protected ABSs are dubbed 
as topological superconductorj^rJ^ and they are charac- 
terized by discrete symmetries such as the particle-hole 
symmetryJ^ii^ Furthermore, it has been clarified as the 
bulk/edge correspondence that when a gapless ABS is 
generated on the edge, the corresponding topological in- 
variant exists in the bulk systemJ^"— 

For quasi-two dimensional superconductors, according 
to the energy dispersions, ABSs are classified into three; 
flat-type, chiral-type and helical type. A fiat-type ABS 
is protected by a one-dimensional winding number which 
is defined for a fixed ky, where ky is the momentum par- 
allel to the surfaceJii^ This fiat type ABS is dubbed 
as mid gap Andreev bound state and is generated for 
nodal superconductors like spin-singlet dxy-wave one or 
spin-triplet p^^-wave one. A chiral-type ABS has a linear 
energy dispersioni^Ti^^ and is realized in spin-triplet chi- 
ral p-wave superconductors like Sr2Ru04i^i^ This ABS 
is protected by the Chern number— i^^i^^i^^ and it causes 
a spontaneous current along the surface. A chiral ABS 
has been recognized as a chiral Majorana edge mode if 
the spin degree of freedom is quenched<^ 

By contrast to the case of chiral ABSs, the time rever- 
sal symmetry is preserved in helical ABSs. A helical ABS 
supports two linear energy dispersions with opposite ve- 
locities, which form a Kramers pair. Instead of a charge 
current, a spin current is spontaneously generated along 
the edge. Stability of the helical ABS is guaranteed by 
the Z2 topological invariant, as in the case of quantum 
spin Hall insulators (QSHs),— , and they are expected to 
be realized in non-centrosymmetric superconductors^^"— 
and a bilayer Rashba system^^^ where the spin-orbit cou- 



pling is important. Several new features of helical ABSs 
have been predicted)^"— Surface ABSs in three dimen- 
sional systems have been also studied. A cone-type ABS 
is predicted for a superconducting analogue of the super- 
fluid ■^He B phase. This ABS is interpreted as Majorana 
fermion^i^"— obeying massless two dimensional Dirac 
equation. Moreover, surface ABSs with various com- 
plicated energy dispersions^ appear in superconducting 
topological insulators?^ e.g., Cua;Bi2Se3^i^. 

A new direction for realization of Majorana fermions 
is to fabricate topological superconductors with conven- 
tional pairingi^ Especially, much attention has been paid 
to realize chiral Majorana modes from the view point 
of topological quantum computing)^"— There are sev- 
eral proposals to fabricate Majorana fermions in systems 
coupled to superconductor via the proximity effect. It 
has been proposed that a chiral Majorana edge mode is 
produced at the interface of ferromagnet / spin-singlet s- 
wave superconductor junction on the substrate of three- 
dimensional topological insulatori^"— Also, a simpler 
scheme using the Rashba spin-orbit interaction and the 
Zeeman field has been proposed.—"— The essential point 
is the simultaneous presence of the strong spin-orbit cou- 
pling and the time reversal symmetry breaking by the 
Zeeman field. There is another way to realize chiral Ma- 
jorana edge modes by using chiral edge states of a quan- 
tum anomalous Hall system (QAH)i^ A QAH can be 
realized by doping of magnetic impurity in a QHSi^ In 
this scheme, the presence of the chiral Majorana edge 
modes can be controlled by the band mass m, chemical 
potential fi and the pair potential A. The number of chi- 
ral Majorana edge modes can be classified by the Chern 
number TV of the systemj^ 

Stimulated by the idea of Qi et al^^ in our previous pa- 
per, we have calculated the energy spectrum of the edge 
states and the resulting surface local density of states 
(SLDOS) for various values of the Chern number Af in 
a heterostructure of a QAH and a spin-singlet s-wave 
superconductor {QAH-\-s)^ To clarify the difference be- 
tween the TV = 1 and TV = 2 states, we applied Zeeman 
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magnetic fields. We have found that when the direction 
of the magnetic field is parallel to the interface, the de- 
generacy of the two chiral Majorana edge modes in A/" = 2 
states is lifted. We have also clarified that the degener- 
acy is lifted by shifting the chemical potential from zero. 
Although the SLDOS has been calculated in detail, the 
relevance to the actual tunneling conductance observed 
in QAH-l-s system have not been clarified yet. 

The purpose of this work is to present a theory of the 
tunneling conductance in this system. If the ABS has a 
flat dispersion, which is realized in high-Tc cuprate, the 
tunneling conductance is expressed by the SLDOS.— 
In the present case, however, the correspondence is not 
clear. Since the ABS has a linear dispersion, the SLDOS 
does not always coincide with the tunneling conductance 
in normal metal (N) / superconductor junction even in the 
low transparent limit.— 1^>2^"— A similar situation occurs 
in three dimensions. Differently from the case of the su- 
perconducting analogue of '^Hc B phased the tunneling 
conductance for the junction of N/superconducting topo- 
logical insulator shows a single zero-bias peak by taking 
into account a finite temperature effect^S or transmis- 
sivity at the interface, even though the SLDOS has a 
double peak structure.— Because of the difficulty to pre- 
dict the charge transport property from the SLDOS, as 
mentioned above, we have to calculate the tunneling con- 
ductance of N/ (QAH-f-s) /N junction by explicitly solving 
the Bogoliubov-de Gennes (BdG) equation. 

The organization of the paper is as follows. In Sec. 
ITU we review the model of QAH with spin-singlet s-wave 
superconductor. In addition, we formulate the tunneling 
conductance in N/ (QAH-|-s)/N junction using the recur- 
sive Green's function. In Sec. IIIIl we calculate the energy 
dispersion, the SLDOS and the tunneling conductances 
in N/QAH/N junction and N/(QAII-hs)/N one. We re- 
veal that an even-odd effect in the angle-resolved con- 
ductance occurs due to interference of Majorana fermions 
whereas the corresponding SLDOS does not. In Sec. IIVI 
we summarize our results. 



II. FORMULATION 

In this section we show the model Hamiltonian of 
QAH+s and the method of numerical calculation for the 
SLDOS. The model of N/(QAH-hs)/N junction and the 
formula of the tunneling conductance with the recursive 
Green's function are also shown. 



A. Hamiltonian of QAH-l-s 

We consider a QAH on the two-dimensional square lat- 
tice, which is obtained by the replacement kx,y — > sin k^^y 
and kl + ky ^ 4- 2(cos 

kx ~t" cos Ajy) in the model used 
in Refs.— 1^ Near the T point, this replacement does not 
change the low energy and low wavelength physics of the 
system. Compared to the continuum model, the square 



lattice model is convenient when we calculate the SL- 
DOS. In the momentum space, the Hamiltonian has the 
form as HqahC^) = d{k) ■ s with 

d{k) = {A sin kx, Asm ky,m{k)), (1) 

where Si is Pauli matrix in spin space and m(fc) = 
TO + 2B{2 — coskx — cosfcy). The band mass term TO(fe) 
determines the magnitude of the energy shift between 
up and down spins. A, B and m are material param- 
eters corresponding to the velocity of the surface Dirac 
fermion, the inverse effective mass of conduction/ valence 
bands, and the band gap, respectively. The sign oim/B 
determines the topological property of the system. Here 
note that the presence of B term is crucial to exhibit a 
QAH. The energy dispersion of the above Hamiltonian is 
symmetric with respect to the mass term m for i? = 0, 
but is asymmetric for B^O. In other words, a nonzero 
value of B makes the sign of to meaningful. Hereafter, 
we take A = B = 1 and the lattice constant being unity 
in our calculations. 

In the following, we consider the proximity effect by an 
attached spin-singlct s-wave superconductor, where the 
pair potential is induced in the QAH (hereinafter we refer 
to it as QAH-|-s). The system is described by the BdG 
Hamiltonian, 

HBdG(fc) = d^{k)s, + [dx{k)sx + dy{k)sy\T^ ~ iiT^ + /\tx, 

(2) 

where Ti is Pauli matrix in Nambu space, /x is the chemi- 
cal potential, and A is the induced pair potential of spin- 
singlet s-wave superconductor. The energy gap of T^BdG 
at fc = is given hy = \m\ — ■^/A^TF/?. The present 
system has three phases, i.e. Af = 0,1, and 2 phases,— 
which arc realized in to > + /.(^ , |?7i| < yA^-T/?, 

and TO < — -^Z A^ + fj? , respectively. 

B. Surface local density of states 

In order to obtain the SLDOS at the edge {x = 1), 
we introduce an infinite potential barrier at a; = 0. We 
calculate the Green's function at a:; = 1 by t-matrix 
formulationi^ The system is infinite along the y-direction 
while it is semi-infinite along the x-direction. Since trans- 
lational invariance is absent along the x-direction, only 
the momentum ky in the y-direction is a good quantum 
number. We express the Green's function in the spatial 
coordinates x and x' for fixed ky as follows; 

Gxx'(ky,ui) = gxx'(ky,Lo) 

- 9xo{ky , uj)goo {ky , uj)gox' {ky,uj), (3) 

with 

gxx'{ky,uj) ^ ^^e'^-^^-'''^g{kx,ky,uj), (4) 

and 

g~'^{kx,ky,uj) ^ uj -'HBdG{kx,ky), (5) 
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FIG. 1. (color online) Schematic illustrations of N/QAH/N 
(a) and N/(QAH+s)/N (b) junctions . 



with 

t{ky) = [m + 2i?(2 — co&ky)\sz + ^sinfeyS^T^ + Ar^. , 

(10) 



and 



^QAH = -Bs^ - 



(11) 



The Hamiltonians of normal electrodes located in the 
left (Hl) and the right (Hr) regions arc given by 



Hhiky) ^ ^ cl{ky){~2tjq cos ky - fj.jq)T^Cx{ky) 

X— — 00 

+ \Y1 4(fcy)(-iNKc,+i(/cj,) + h.c. ) , (12) 



where is the total number of lattice points in the x- 
dircction. In the right hand side of Eq.(|31), the first term 
denotes the unperturbed bulk Green function, and the 
second term comes from the scattering effect at the edge. 
The angle resolved SLDOS N{ky,uj) at x = 1 is written 
as 



where 



N{ky,uj) = --ImTr [P,Gfi{ky,uj)] , (6) 
G^,,{ky,Lo) = Gr.'iky,^ + ir,), (7) 



is the retarded Green's function, r/ is an infinitesimal pos- 
itive number, and Po = (l + ''"2)/2 is the projection opera- 
tor onto the particle subspacc. From the above equations, 
one obtains the SLDOS D{lo) for energy tu measured from 
the Fermi level as follows 



(8) 



where Ny is a total number of lattice points for the y- 
direction. In the actual calculation, we set Nr^ ~ Ny ~ 
4096. 



C. N/(QAH+s)/N junction 

Now we show the Hamiltonian of N/ (QAH-|-s) /N junc- 
tion illustrated in FigH] The Hamiltonian of QAH+s in 
the center region of the junction is represented as 



Nc 



x=l 

A^QAH-l 

+ H{ky)tQAliC, + i{ky)+h.C.], (9) 



'Hniky) ^ ^ cl{ky){-2t-!^ cos ky - fiT^)T^Cx{ky) 
oo \ 

cl{ky){-tN)Tza^+i{ky) +h.c.\ . (13) 

yX = NQAll + l / 

We also assume the following simple hopping "Hj between 
the electrode and the QAH+s, 

'^iiK) = 4(fcy)(-^j)^^ci(A:y) 

+ cjvQAjj(^-y)(-ij)T2CjVQAH+l(fci;) + h-c- (14) 

In the actual calculation, tj is fixed as tj = tN, for sim- 
plicity. 



D. Tunneling conductance and 
recursive Green's function 



The angle resolved tunneling conductance ^{py,u]) in 
a junction is given by the Lee-Fisher formulap^^ 

~ ^X,X^X+l,X+l ^ ^X+l,X+l^X,x)]7 (15) 

with G"^., — Im G^^, . Due to current conservation in the 
normal metals, we can choose arbitrary x for Eq. (fT5l) in 
X < or a; > Nqah + 1, except in the superconducting 
region {1 < x < iVqAn)- The total conductance T is given 
by r(w) = J2p liPyj^)- We first calculate the Green's 
function Gh,i,i in the left semi-infinite system where the 
sites in a; > / are deleted. Gi,,x.x satisfies the following 
recursive relation: 

Gl,!-,^; ~ 9x^ ^ 'Hx.x~iGh,x-i,x-i'Hx-i,x, (16) 

with 'Hx,x' being the hopping from x' to x. Here, 
g~^(jpy,ui) — Lo ~ T~Lx.x{Py) is the Green's function in the 
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isolated x-th column. In the present model, only "Hx.x' 
with |a; — x'l < 1 is nonzero, and given by 

'Hx,x' = [(-2tN cos ky - I-Ln)Sx,x' ~ tT<!5x',x±l]Tz, (17) 

for a;, a;' < or x, x' > A^qah + 1, and 

y-x.x' = <^{ky)5x,x' + ^QAH'^2;',2; + l + t^Q^p^Yl^x' ,x-l, (18) 

for 1 < x,x' < Nqah- At the interfaces, 'Ho.i, Hi^, 
and Hatqah+i.Wqah arc given by 



dispersion 



SLDOS 



conductance 



n 



t 



H + 1 — ^JVqah + 1,Wqah 



(19) 



The Green's fmiction in the right semi-infinite system 
Gb.,1,1 where the sites in x < ^ are deleted satisfies the 
following relation. 



^K.x X " 9x^ ^ T~l-x,x+lG^^x + l,x + lT~i-x+l,. 



(20) 



It is noted that the Green's functions at the edge of the 
electrode (Gl,o,o and Gr,^Nqah+i,Nqah+i) are obtained 
by Eq.jSl). Then, using Eqs. ^ and we can re- 

cursively obtain Gi^^x,x and G^^x,x for any x. The site- 
diagonal part of the Green's function is obtained in terms 
of the above Green's functions, 

Gx^x ~ 9x^ ^ T~Lx,x-iG\^^x-i,x-iT~Lx-i,x 

^ T~l-x,x+iG-R,^x+i,x+iT~ix+i,x, (21) 

and the site-off-diagonal parts are also obtained as 

Gx,x+1 — Gx,xT~(-x,x + lGT;i,x+l,x + l7 (22) 

Gx+i,x = Gx+i.x+i'Hx+i,xGi^^x,x- (23) 

We can calculate the conductance T from Gx,xi Gx+i,x+i, 
Gx.x+i, and Gx+i,i by using Eq. ([151). 



III. RESULTS AND DISCUSSIONS 

In this section, we show our numerical results for the 
electronic states and the tunneling conductance of the 
QAH-l-s. Experimental proposals to detect our results 
are also discussed. 



A. Electronic states and tunneling conductance in 
N/QAH/N junction 

Before discussing the superconducting case, we check 
the electronic states and the tunneling conductance of the 
N/QAH/N junction at A = ^ = 0. The energy disper- 
sions with the finite width (A*qah = 100) are shown in 
Figs. [2ja), (b), and (c). There is a band gap in Figs, 
and[2ljc), respectively, while it closes at the critical point 
with m~Q [Fig. djb)]. Gapless chiral edge modes show 
up when m < [Fig. HJc)]. Note that two gapless modes 
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FIG. 2. (color online) Energy dispersion relations [(a),(b), 
and (c)], SLDOS [(d), (e), and (f)], and the tunneling con- 
ductances in N/QAH/N junction, for the case with the trivial 
insulator [{g):m/B — 0.5], the critical point [(h):m — 0], and 
the QAH [{c):m/B — —0.5]. a denotes the lattice constant. 
The parameters are taken as follows. A/B = 1, A = n = 0, 
and A^QAH ~ 100 for the energy dispersion, A^qah ~ 4096 for 
the SLDOS, A'qah = 1000 for the conductance. 



propagating in opposite directions appear in Fig. HJc) 
since both the left {x = 0) and right edges {x ~ A^qah) 
are present in the calculation, i.e., each edge has an edge 
state. Figures Eljd) , (e), and (f) show the SLDOS at the 
edge. In the trivial insulator phase with m > [Fig. 
^d) and (e)], the line shapes of the SLDOS are the same 
as those in the bulk. On the other hand, in the case of 
QAH, the SLDOS is enhanced in the band gap due to 
the gapless edge modes, as shown in Fig. [IJf). 

The line shapes of tunneling conductance are similar 
to those of the SLDOS when m > [FiglS^g) and (h)], 
i.e.; U-shaped gap in |a;| < m [Figl^Jg)] and V-shaped dip 
at m = [Figlljh)]. In the case of QAH, ahhough the 
SLDOS shows a zero-energy peak [FiglSl^f)]. the corre- 
sponding tunneling conductance shows a U-shaped gap 
[FiglDJi)] similar to that in the case with to > [Fig. 
[2l^g)]. This is because the central region of the junction 
has a bulk gap, thus the tunneling conductance should 
be zero. 

Here we notice that the conductance in Figs. [Hg) and 
(i) takes a small but non-zero value in the energy gap, 
but this comes from a nonzero value of 77 in our numerical 
calculation. Indeed, as one decreases 77 in Eq. (O and 
increases the number of QAH layer (TVqah), 7(0, 0) tends 
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Nqah n 

FIG. 3. (color online) Angle resolved tunneling conductance 
')(ky, oj) at ky = and cj = as functions of A'^qah (a) and 77 
(b), for the QAH with m/B — —0.5 and the trivial insulator 
with m/B = 0.5 and A^^qah = 1000. The parameters are the 
same as in Fig. [21 

to be zero, as shown in Fig. [31 i.e., the corresponding 
conductance goes to zero in the limit of A^qah oo 
and rj — > +0. Additionally, we note that the tunneling 
conductance of the QAH is lager than that of the trivial 
insulator for finite iVqAH and ij (Fig [31). This is due to 
hybridization of the edge states located at x' = and 
X = iVqAH- 

B. Electronic states of QAH+s 

Now we consider the superconducting case with 
A/B ~ 0.25. The trivial superconductor with TV = 
does not have any gapless state as shown in Fig. 
mja), while the topological superconductors with TV = 1 
[Figgi;b)] and TV = 2 [Figgljc)] have. Although the gap- 
less modes in the latter two phases have a similar energy 
dispersion, we can distinguish them by the SLDOS, as 
shown in Figs. [HJe) and (f): The line shape of the SL- 
DOS for TV = 1 shows a zero-bias peak [FigHl^e)]. On 
the other hand, that for TV = 2 shows a larger zero-bias 
peak and satellite peaks at oj/Eg ^ ±2 [Fig. [HJf)]. The 
larger zero-bias peak is due to two gapless modes, and 
the satellite peaks come from a branch of ABS near the 
bulk bands. 



C. Even-odd effect in the angle resolved 
conductance 

Next we discuss the angle resolved conductance 
j{ky,uj) in the N/(QAH+s)/N junction shown in Fig. 
[3 In TV = [Fig. ISJa)] and TV = 1 [Fig. [S^b)] phases, 
the conductance spectra are naturally understood by the 
energy spectra [Figs. IHa)(b)] and the SLDOSs [Figs. 
[3l^d)(e)]. Due to resonance between the incident state 
and the chiral edge mode, the value of conductance in 
TV = 1 phase takes ^{ky,uj) ^ 2e^/h at cj = -Eodgo(fcy), 
where i?odgo(fcy) is the energy dispersion relation of the 
edge state. Note that the conductance spectra are asym- 
metric with respect to fcj, because the present edge 
mode is chiral. 




-1 1 0.05 0.1 

ky [q-1] SLDOS 

FIG. 4. (color online) Energy dispersion relations and the 
SLDOSs in QAH+s with M = Q [(a) and (d): m/B = 0.5], 
TV = 1 [(b) and (e): m = 0], and TV = 2 [(c) and (f): m/B = 
—0.5]. The pair potential is taken to be A/B = 0.25, and the 
other parameters are the same as in Fig. [21 



Angle resolved conductance 'y{ky,u:) 
(a) TV = (b) .V = 1 (c) A' = 2 




-1 -1 -1 1 

ky [a~^] ky [a^^] ky [a^^] 



FIG. 5. (color online) Angle resolved conductance in 
N/(QAH-|-s)/N junction with TV = (a), TV = 1 (b), and 
TV = 2 (c) phases. 



In TV = 2 phase, we obtain a remarkable result: In this 
phase, it is natively expected that the tunneling conduc- 
tance take a doubled value of that in TV = 1 phase since 
there are two edge modes. It is, however, not the case. 
As shown in Fig. [Sjc), in TV = 2 phase, the conduc- 
tance takes a smaller value than that in TV 1 phase at 
u! ^ Ecdgc{ky)- In particular, the conductance vanishes 
just aX OJ ~ Ecdgc{ky)- To confirm this, we study j{ky,iu) 
at ky = and w = as functions of TVqah [Fig- El^a)] 
and 77 [Fig. Wlb)] . It is found that as TVqah increases and 
rj decreases, 7(0,0) converges to 2e^/h in TV = 1 phase. 
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Wqah t] 

FIG. 6. (color online) Angle resolved conductance 7 at fcy = 
and (xJ = as functions of A'^qah (a) and r\ (b). 




FIG. 7. (color online) The energy dispersion of the QAH+s 
(a) and the angle resolved conductances (b) in A/" = 2 phase. 
The chemical potentials are set to be At/_B = 0.3. The energy 
gap is given by E^/B = 0.11. 

and to in A/" = 0, 2 phases. Therefore, the conductance 
shows an even-odd effect as 7 = [1 — (— l)^]e^//i in the 
presence of chiral Majorana fermions. 

The vanishment of the tunneling conductance origi- 
nates from the degeneracy of Majorana edge fermions in 
Af = 2 phase: Indeed, it is suppressed by lifting the de- 
generacy by tuning the chemical potential or applying 
Zeeman fields. At a finite chemical potential, the degen- 
eracy of the gapless modes is lifted and the zero energy 
states appear at the two finite ky points, as shown in 
Fig. [Tl^a). The conductance at w = with the finite /x is 
shown in Fig. Eljb). Contrary to that with /U = 0, 7(0, 0) 
remains non-zero. This result suggests that the vanished 
conductance arises from interference of degenerated Ma- 
jorana fermions. 

Here we note that the degeneracy of Majorana edge 
fermions is ensured by a symmetry of the system: When 
/J, = 0, the BdG Hamiltonian ^ has the following addi- 
tional symmetry, 

SyTzHBdcikx, ky)syTz = -HBdcikx, -ky), (24) 

and on the axis, this reduces to the so called chiral 
symmetry. 

{r,HBdG(fc.,0)} = 0, (25) 

with r = SyTz- Thus, following RefsllTl andlssl one can 
introduce the one-dimensional winding number, 

^^"i //fc.tr[rHB^c5,,?^BdG],^^o' (26) 



which can be evaluated as = 2 in the case of A^ = 
2.—. Therefore, the bulk-edge correspondence ensures 
that there exist two degenerate Majorana edge modes 
at ky = 0. As we mentioned above, because the van- 
ishment of the tunneling conductance occurs only when 
Majorana edge modes are degenerate, it is very likely 
that this chiral symmetry is responsible for the destruc- 
tive interference of the tunneling conductance reported 
here. 

In the normal (not superconducting) states, even- 
odd effects in conductance appear in mono/bi-layer 
graphene?^ which can be generalized to the system 
with spin-orbit interactionsi^ They are interpreted as 
a result of mirror symmetry of the systemj^ Also, 
graphene nano-ribbons show even-odd effects in the 
conductances^"— , which can be understood using parity 
of the system. We believe that our result is also explained 
in the viewpoint of symmetry. 

On the other hand, in superconducting states, vari- 
ous even-odd effects of Majorana fermions have been re- 
ported so far. For instance, in N / a chain of Majo- 
rana bound states junction, the tunneling conductance 
shows an even-odd effect as a function of the length of 
the chaiui^ It has been also known that the SLDOS at 
the zero energy in a multiband Rashba superconductor 
with Zeeman interaction shows an even-odd effect as a 
function of the number of occupied subbandi^ In these 
cases, no degeneracy of Majorana fermion exists by hy- 
bridization when the number of the Majorana fermions is 
even. On the other hand, in our case, the degeneracy in 
TV = 2 phase is essential to obtain the even-odd effects. 
Therefore, the even-odd effect reported in the present pa- 
per is essentially distinct from the previous ones, and it 
originates from the interference without using interfer- 
ometers as discussed in Refs. [H, [s^ andl92l 



D. Proposals for experiment 

Before closing the section, we propose how to detect 
the even-odd effect mentioned above. The simplest ob- 
servable is the angle-integrated tunneling conductance. 
In TV = phase, the system has no gapless mode, then 
the value of conductance becomes smaller in the energy 
gap, as shown in Fig. [8] (red line). In TV = 1 phase, the 
line shape of conductance shows a zero-bias peak due to 
gapless mode [Fig. [5] (green line)]. On the other hand, 
in TV = 2 phase, the line shape of conductance shows a 
zero-bias dip [Fig. |8](blue hue)], in spite of the presence 
of gapless modes. 

The more direct evidence is to observe the angle re- 
solved current by scanning a charged tip above the 
system as was performed for two-dimensional electron 
gases in GaAs heterostructures.— Moreover, it is use- 
ful to fabricate the nanoribbon of N/(QAII-|-s)/N junc- 
tion since normal incident electron with ky = effectively 
contributes to 7(0,0). 
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FIG. 8. (color online) Tunneling conductances in the 
N/(QAH+s)/N junction with A/" = 0, A/" = 1, and A/" = 2 
phases. 



IV. SUMMARY 

We studied the tunneling conductances of 
N/(QAH+s)/N junction in A/" = 0, A/" = 1, and 
N ~ 2 phases. In the presence of gapless edge modes in 
N ~ 1 and N = 2 phases, the corresponding SLDOSs 
take finite values in the bulk energy gap. In A/" = 1 
phase the angle resolved conductance also takes the finite 
value of ■~f{ky,uj) = 2e^/h when the incident electron is 
resonant with the chiral edge mode at uj = Ecdgc{ky)- On 
the contrary, in A/" = 2 phase, the tunneling conductance 
vanishes at w = i?edgo(^y) although the corresponding 
SLDOS does not. This stems from the interference of the 
degenerated Majorana fermions. Namely, an even-odd 
effect with respect to the number of Majorana fcrmion 
N occurs. 

Although we have partly addressed the mechanism of 
the even-odd effect, it has been not definitely answered. 
To reveal this, the following things are needed to be un- 
veiled: relation between the chiral symmetry and the 
tunneling conductance, robustness of the even-odd effect 
against disorder proved by the microscopic calculation, 
and the even-odd effect for the higher Chern number of 
A/" > 3. Wc will study these issues in the future work. 



Let us consider a normal metal in the left side < 0), 
whose Hamiltonian is given by 



2mr. 



(Al) 



and QAH-l-s in the right side {x > 0). Here k = 
ik'^ + kyY^^ is the magnitude of the two-dimensional mo- 
mentum. The Hamiltonian of QAH-l-s is obtained by the 
k ■ p theory as^S 

H{k) = m{k)sz + Aa{kxSx + kySy)Tz - + Ar^, 

(A2) 

m{k) = m. + B{kaf . (A3) 
The eigenvalue of the above Hamiltonian is given by 

Eo,p{k) = a m^{k) + A^k^a^ + + 

+ 2(3^J{m?[k) + A^k^a^)n^ + m2(fe)A2j , (A4) 

where a,/3 = ±. The corresponding eigenvector Uap{k) 
is also obtained analytically. 

Now we calculate the tunneling conductance, gen- 
eralizing theories of the tunneling spectroscopy of 
conventional^^ and unconventional^i^ superconductors. 
The wave function in the normal metal {x < 0) is given 
by 



i>N,six) 



+ ass'Xs'he 



Jkyy 



(A5) 



where Xst is the eigenvector of H^[k) with spin s 
for electron (r = e) or hole (t = h), and k^x = 
{kl-kl^/'^ = fee cos 61, ke = ^/2rnJJi^r+E), khx = 
[2mo(/XN ~ E) — kyY^'^. The first term of the wave func- 
tion denotes an injected electron, and the second (third) 
one denotes a reflected hole (electron) with reflection co- 
efficient a^gi {hggi). The wave function in the QAH-|-s 
{x > 0) is given by 



i{qiX+kyy) 



(A6) 
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Appendix A: Tunneling conductance of N/(QAH-|-s) 
junction in the continuum limit 

To confirm the even-odd effect found in this paper, we 
study the tunneling conductance of N/ (QAH+s) junction 
in the continuum limit. 



where g^, (i = 1, • • • , 4) is a solution of i? = Eo^iSi [qi, ky). 
Among the eigenvectors, VQAH-i-s(a') consists of those 
with Eaii3i{qi,ky)/dqi > or hn{qi) > 0, where 
the former denotes right-going states and the latter 
denotes localized states in the vicinity of a; = 0. 
These wave functions are connected at the interface 
{x = 0) by the conditions^l, Vn(0) = ?/'qah+s(0) and 
VNtp^iO) = i'QAH+sV'QAH+s(0), with the velocity opera- 
tor 'i;n(qah-hs) = ^H^|^QAR+s)/^kx\k^-^-i^^■ The above 
equations determine the coefficients Uss', bss' and U. Fi- 
nally, the charge conductance ^{ky,bj) is given by 



l{ky,u) 



(A7) 
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infinitesimal small factor ry is not necessary. The above 
2 result indicates that the even-odd effect found in this 
paper is robust, and the vanishing conductance is driven 
^■■^ by the degenerating two Majorana fermions. 

1.0 

0..5 


-10-10 1 
ky [a^-'] ky [a^^] 



FIG. 9. (color online) Angle resolved conductance 7(fcy,tj) 
of the N/(QAH+s) junction for A/" = 1 with fl = fi/B = 
(a) and jv = 2 with fl = 0.3 (b) phases. E^, = 0.25-B is the 
magnitude of the superconducting gap. 




FIG. 10. (color online) The angle resolved tunneling conduc- 
tance y{ky, 0) at the zero-bias voltage uj — for J\f = 1 and 
Af = 2 phases (a). Those in M = 2 phase for p. = fi/B = 
and jl — 0.3 are also shown (b). 



In the following, the material parameters of the normal 
metal are fixed as rricBa? = 1, fin/B = 100, and the 
material parameters of QAH-l-s are chosen a.s A = B 
and A/B = 0.25. 

The obtained angle resolved tunneling conductances 
are shown in Fig. [5] These spectra are consistent with 
those obtained in the lattice model shown in Fig. [S] the 
value of 7(/cy, £^edge(fcy)) takes 2e^//i in A/" = 1 phase 
[Fig. mja)] while it takes in TV = 2 phase [Fig. Mh)]. 
In order to see the even-odd effect more clearly, we fo- 
cus on 7(fcy,0). Figure fTUTa) shows the angle resolved 
conductance at the zero-bias voltage (w = 0), where the 
branch of Majorana fermions appears at ky = 0, i.e., 
E^dgeiO) = 0, both for A/" = 1 and A/" = 2 phases. The 
value of 7(0,0) takes 2e^/h in A/" = 1 phase and in 
Af = 2 phase. However, if one tunes the chemical po- 
tential fl away from zero, where the degeneracy of two 
Majorana fermions is lifted, the value of conductance re- 
covers to be finite at ky = 0, as shown in Fig. [TUTb). 

As compared to the calculation in the lattice system, 
the present approach in the appendix has advantages, 
i.e., it is easy to take the thermodynamic limit, and the 
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